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Abstract 

in 

! We provide detailed arguments on how to derive properties of generalized form 

I factors, originally proposed by one of the authors (M.K.) and Weisz twenty years 

ago, solely based on the assumption of "maximal analyticity" and the validity of the 
LSZ reduction formalism. These properties constitute consistency equations which 
CIhI allow the explicit evaluation of the n-particle form factors once the scattering matrix 

^ ' is known. The equations give rise to a matrix Riemann-Hilbert problem. Exploiting 

the "off-shell" Bethe ansatz we propose a general formula for form factors for an 
odd number of particles. For the Sine-Gordon model alias the massive Thirring 
^ ' model we exemplify the general solution for several operators. In particular we 

. calculate the three particle form factor of the soliton field, carry out a consistency 

check against the Thirring model perturbation theory and thus confirm the general 
formalism. 
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1 Introduction 



More than fifty years ago, Heisenberg pointed out tlie importance of studying analytic 
continuations of scattering amplitudes into the complex momentum plane. The first 
concrete investigations in this direction were carried out by Jost and Bargmann 0], 
initially for non-relativistic scattering processes. The original ideas turned out to be very 
fruitful and lead to interesting results on-shell, i.e. for the S-matrix |Q, as well as off-shell, 
that is for the two-particle form factors, see for instance ^ . 

Once one restricts ones attention to 1+1 dimensional integrable theories, the n-particle 
scattering matrix factories into two particle S-matrices and the approach, now usually 
referred to as the bootstrap program, reveals its full strength. On-shell, it leads to the 
exact determination of the scattering matrix |^, (for reviews see also [9-12]). The 
results obtained in this way agree with the S-matrix obtained from the extrapolation of 
semi-classical expressions for the Sine-Gordon model [Q. The first off-shell considerations 
were carried out about two decades ago by one of the authors (M.K.) et al. |TB|, 0, who 



introduced the concept of a generalized form factor and formulated several consistency 
equations which are expected to be satisfied by these objects. Thereafter this approach 
was mainly developed further and studied in the context of several explicit models by 
Smirnov et al. [15-23]. Recently this program has seen some revival in relation to models 
which arise as perturbations of certain conformal field theories particularly in the 
context of affine Toda theories [^ and closely related models [26-48] . 

An entirely different method, the Bethe ansatz was initially formulated in order 
to solve the eigenvalue problem for certain integrable Hamiltonians. The approach has 
found applications in the context of numerous models and has led to a detailed study of 
various mass spectra and S-matrices (for reviews and an extensive list of references see 
for instance [0]). The original techniques have been refined into several directions, of 
which in particular the so-called "off-shell" Bethe ansatz, which was originally formulated 



by one of the authors (H.B.) [^T], will be exploited for our purposes. This version 
of the Bethe ansatz paves the way to extend the approach to the off-shell physics and 
opens up the intriguing possibility to merge the two methods, that is the form factor 
approach and the Bethe ansatz. The basis for this opportunity lies in the observation 
)3, |55[ , that the "off-shell" Bethe ansatz captures the vectorial structure of Watson's 



equations (see section ^]2| properties (i) and (ii)). These are matrix difference equations 
giving rise to a matrix Riemann-Hilbert problem which is solved by an "off-shell" Bethe 
ansatz. Furthermore, there exist interesting speculations in order to make contact with 
general concepts of algebraic quantum field theory [6^, |61[] . 

Conceptionally, the on and off-shell approaches are very similar. For the on-shell 
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situation one has certain constraints resulting from general physical and in particular 
analytic properties (referred to as "maximal analyticity assumption"), which lead to a 
set of conditions which turned out to be so restrictive that they allow to construct the 
exact scattering matrix almost uniquely. This approach is adopted in order to determine 
the key off-shell quantities, i.e. the form factors. In the present manuscript we provide a 
detailed derivation of the consistency equations solely based on the maximal analyticity 
assumption and the validity of the LSZ-formalism |^ (see also Form factors are 



vector valued functions, representing matrix elements of some local operator 0{x) at the 
origin between an in-state and the vacuum, which we denote by (refer equation (|3.2|) for 
more details) 



[iiP^ + PjY + ^eh<^<J<n)) :=(0|O(0)|pi,...,pX....a„ • (1-1) 

Once all the ra-particle form factors are known, one is in principle in a position to compute 
all correlation functions. In particular the two point function for an hermitian operator 
O in real Euclidian space reads 



{Oix)Om = t f^^\j-^\F^ie,,...,e^) pexp (-r±m,coshe- 



1.2) 



Here r denotes the radial distance r = \Jxi + and 9 is the rapidity related to the 
momentum via Pi = m^sinh^j (see section 3.2 for more details). The explicit evaluation 
of all integrals and sums remains an open challenge for almost all theories, except the 
Ising modelQ. Important progress towards a solution of this problem has recently been 
achieved in []59| . 

A commonly used procedure which will yield expressions which satisfy all of the con- 
sistency requirements is constituted out of the following steps: First of all one has to 
have solved the on-shell system, that is one requires expressions for the S-matrix. In the 
next step one usually makes an ansatz for the form factors of a type already introduced 
in 



)13| , in which one extracts explicitly the expected singularity structure. The nature 
of the ansatz guarantees by construction that the generalized Watson's equations (prop- 
erties (i) and (ii)) are satisfied once the scattering matrix is diagonal. For generically 
non-diagonal scattering matrices one may invoke also the techniques of the "off-shell" 



Bethe ansatz |53, 54, 55 1 in order to capture the vectorial structure of the form factors. 



^Of course one may also adopt a very practical point of view and resort to the well-known fact 



that the series expansion of correlation functions in terms of form factors (L2) converges very rapidly. 
Consequently correlations functions may be approximated very often quite well by simply including 
the two-particle form factor into the expansion. /From that point of view the form factor program is 
completed, since the calculation of the two-particle form factors is well understood. 
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The ansatz only involves the rapidity differences, apart from a possible pre-factor, which 
takes the spin of the local field O into account, and has therefore the desired behavior 
under Lorentz transformations (refer property (v) in section 3.2). General solutions for 
the so-called minimal form factors (the function which satisfies the functional equation 
(4.10)) are always fairly easy to find. Once the scattering matrix is non-diagonal one has 
also to encode the vectorial structure at this stage. What is then left, is to determine 
a general function which takes the complete singularity structure into account. For this 
purpose one may now invoke properties (iii) and (iv) (equations (3.12) and (3.13) for the 
bosonic case), which lead to a set of recursive equations. In principle these equations 
may now be solved step by step, once the first non-vanishing form factor for a particular 
operator is properly fixed. However, only after a few steps the expressions become usually 
algebraically very complex and reveal very little insight. Therefore, it is highly desirable 
to search for structures of a more general nature, that is in particular to seek for closed 
expressions for all n-particle form factors. Only such expressions may ultimately shed 
more light on the analytic expressions for the correlation functions ( |1.2|) . Alternatively, 
one may try to construct directly a representation for the creation operators of the parti- 
cles in the in-state in (|1.1|) |[6^ , |63t |6^ , |65|| . Representing the local operator O in the same 
space, one may in principle also compute the form factors. 

In the present manuscript we provide a general expression (see theorem 4.1) of a differ- 
ent kind, which solves all the consistency requirements. It is very generic by construction 
and, roughly speaking, captures the vectorial nature of the form factors by means of 
"off-shell" Bethe ansatz states and the pole structure by particular contour integrals. We 
exemplify this general expression for the form factors of the Sine-Gordon model involving 
an odd number of states, which was hitherto unknown. For the even case similar expres- 
sions may be found in [|1^, |1^]. We present a detailed analysis of the three particle form 
factor. 

Once solutions for the set of consistency equations are found, it is highly desirable 
to verify the solutions with some alternative method. Several different methods have 
been developed in recent years. Assuming that the theory under consideration results 
from the perturbation of some conformal field theory, one may carry out the following 
consistency checks. For instance one may take the operator in the form factor to be the 
trace of the energy momentum tensor and exploit the so-called c-theorem |^| in order to 
obtain a first indication about the result. This check is not extremely restrictive what the 
higher n-particle form factors concerns, since the expected value for c is usually already 
saturated after the two-particle contribution. Alternatively one may also compare with the 
perturbation theory around the conformal field theory, which is possible for all operators 
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of the model. The latter approach has turned out to be very fruitful [^. A further 
consistency check consist out of the comparison between the exact result obtained from 
the form factors with the predictions of the renormalisation group (that is asymptotic 
freedom etc. [0, 0). In the present manuscript we present a check of our solutions 
against conventional perturbation theory in standard quantum field theory. 

The manuscript is organized as follows: In section 2 we review the properties of the 
general scattering matrix and in particular the Sine-Gordon S-matrix. In section 3 we 
motivate the general properties of the generalized form factors, for simplicity initially only 
for the bosonic case, which we thereafter extend to the general situation involving also 
fermions. In section 4 we briefly explain the "off-shell" Bethe ansatz and state theorem 
4.1, the main result of the manuscript. We present a general formula^ (based on the 
"off-shell" Bethe ansatz) for form factors with an odd number of solitons or anti-solitons. 
Furthermore, we provide an explicit analysis of several two- and three particle form factors 
and carry out various consistency checks relating different form factors to each other. In 
section 5 we compare our solution for a three particle form factor against perturbative 
perturbation theory. Our conclusions are stated in section 6. In appendix A we provide 
the proofs of the properties of the generalized form factors. In appendix B we proof 
theorem 4.1 and appendix C serves as a depot for several useful formulae employed in the 
working. 

2 The S'-matrix 
2.1 General Properties 

In this section we briefly review some of the well known facts on the general properties of 
the scattering matrices. The Fock space is spanned by the in- or out-states of the particles 

I PU...,Pn = ■ ■ ■ ^t'^'^Pn) I ) (2.1) 

where the a^''s are creation operators. The p's denote the momenta and the a's the 
internal quantum numbers of the particles, such as the particle type etc. We choose the 
normalization 

o,'{p'\p)a = 2uJ 277 S{p' - p) = 477 5(9' - 6) (2.2) 

where the rapidity 9 is related to the momentum by p = m sinh 9 and uj = \/rn? -|- p^. 

In an integrable quantum field theory in l-|-l-dimensions there exists an infinite set 
of conservation laws. Therefore in a scattering process the sets of incoming and outgoing 

^Our formula is similar to an analogous one of Smirnov [jisf for even number of particles. This should 
be a starting point for a comparison of both formulae. 
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momenta are equal 

{Pu---,Pn} = {pi,---,Pn'}- 

The n-particle S-matrix is defined by 

\Pl,---,Pn)Z...a„ = IPl.---.Pn)°t<^^"t£bl,---,Pn) 

The statistics of the particles has been taken into account by the diagonal matrix cr*-"^. It 
is a product of all two particle matrices a with entries —1 if both particles are fermions 
and +1 otherwise (see |^). As a consequence of integrability , i.e. the existence of an 



infinite number of conserved quantities, the n-particle S-matrix factorizes into n{n — l)/2 
two-particle ones 

a(")5(")(pi,...,p„) = n^5(2nP.,P.), 

i<j 

where the product on the right hand side has to be taken in a specific order (see e.g. [0). 
For this reason it is sufficient to investigate the properties of the two-particle scattering 
matrix. As is usual in integrable quantum field theories in 1+1- dimensions it is most 
convenient to regard the two-particle S-matrix as a function of the rapidity differences 
6 = \6i — 6j\ rather than as a function of the Mandelstam variables Sij = {pi + pj)"^- In 
order to establish the analytic properties of the two-particle S-matrix one may employ the 
relations Sij = + + 2mimj cosh 9ij, Uj = {pi —PjY = 2m^ + 2m| — sij. Considering 
the scattering matrix as a function in the complex Sj^-plane, there will be two branch cuts 
present, the s-channel one for Sij > {nii + mjY and the t-channel one for Sij < {m.i — irijY. 
In figure 1 the physical s-channel and t-channel regions in the complex (a) s- and (6) 
6'-plains are labeled by I and II, respectively. The crossing transition is depicted by an 
arrow. This and the transitions corresponding to the exchange of in- and out-going waves 
are given by: 



I ^11 

I ^ in 

II ^ VI 



Sij + ie tij — ie -v^ 9 ^ in — 9 

Sij + ie ^ Sij — ie 9^—9 

tij — ie ^ tij + ie <v4> in — 9 in + 9 



(It will be important in the following to notice that the t-channel cut (II-IV) is not present 
for form factors as a function defined in the complex Sjj-plane.) 

Let ^ be a finite dimensional vector space, whose basis vectors label all types of 
particles of the model. Then one considers the S-Matrix as an intertwining operator 
acting on the tensor product of two of these spaces 

Si2{9) : Vi(E)V2-^V2®Vi. 
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Figure 1: The analyticity domains in the complex planes of (a) the Mandelstam variable 
Sij — ipi o,nd (6) the rapidity difference variable 9 = \9i — 92\. The physical regimes 

in the s- and t-channels are denoted by I and II, respectively. The crossing transition 
from the s- to the t-channel is indicated by the arrow. As explained in the main text, the 
interchange of in and out means transition from I to III for the s-channel, and II to IV 
for the t-channel. The dots denote the possible positions of poles corresponding to one 
particle intermediate states. 



The unitarity of the S-matrix reads 

E (sfS'iO))* S^J{9) = 5^.Jp.p or S,,(-9) S,,(9) = 1 (2.3) 

a'(3' 

since by analytic continuation from positive to negative variable one has Sl2{9) — S2i{—9). 
The crossing relations are 

St}{e) = Sfjtn -9) = Sf,{i7^ - 9) (2.4) 

where the bar refers to the anti-particles. The Yang-Baxter equation which follows from 
the higher conservation laws is 

(<T^)l2(|^12|)((T^)l3(|^13|)(^^)23(|^23|) = (^^)23(|^23|)(^^)l3(|^13|)(^^)l2(|^12|) (2.5) 

where % — 9i — 9j. When there are no transitions of the sort that two bosons change 
into two fermions, the signs given by the statistics cancel. 

As usual we use here and in the following the notation for a vector with components 
yai...an matrix with elements acting on these vector 

e V,...^ = V,^---^Vn, : V,...^ ^ K...n (2.6) 

where all vector spaces Vi are isomorphic to V and whose basis vectors label all kinds of 
particles. An S-matrix as Sij acts nontrivial only on the factors Vi <S>Vj and in addition 
exchanges these factors. If we want to express the fact that a particle belongs to a multiplet 
of a specific type of particles, we also write v"' e Va,v^ e Vj, etc. and consider V — 0„ T4 
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as the direct sum of all these spaces. Usually these spaces Va are the representation spaces 
of a symmetry group or quantum group of the model. 

The physical S-matrix in the formulas above is given for positive values of the rapidity 
parameter 6. For later convenience we will also consider an auxiliary matrix S regarded as 
a function depending on the individual rapidities of both particles 9i, 62 or 6*12 = 61 — O2 



{Sa),,\\9,-92\) for 9, < 92 



(2.7) 



with cr taking into account the statistics of the particles. Up to these statistics factors S 
is obviously the analytic extension of the physical S-matrix S from positive to negative 
values of 9, due to the unitarity (f^.3|). 

It appears convenient to introduce a graphical representation for several of the ampli- 
tudes, which will allow us to develop a more direct graphical intuition for the derivation 
of several relations. The auxiliary matrix S may be depicted as 



S 



121(71 



92) 




Here and in the following we associate a rapidity variable 0j G C to each space Vi which is 
graphically represented by a line labeled by 9i or simply by i. In terms of the components 
of the S-matrix we have g 



9i 

a 




7 

92 

/3. 



In terms of the auxiliary S-matrix the Yang-Baxter equation has the general form 

Sl2{9i2) Si3{9is) 323(923) = >S'23(^23) 'S'i3(6'i3) 5*12(^^12) 

which grphically simply reads 





(2.8) 



Unitarity and crossing may be written and depicted as 



'S'2i(6'2i)5' 



121(^12 j 



1 : 



1 2 



1 2 
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S'i2(^l — 02 



til . 






12 12 12 

where C^^ and Cii are charge conjugation operators with components C"^ = C^^ = 5ap. 
We have introduced the graphical rule, that a line changing the "time direction" also 
interchanges particles and antiparticles and changes the rapidity as 6* ^ iivr, as follows 



— 



e\ \_e-m 

a (3 



a/3 



a P 

e\ le + iw 



(2.9) 



Similar crossing relations will be used below to investigate the properties of form factors. 
2.2 Bound states 

Let the two particles labeled by 1 and 2 of mass mi and m2, respectively form a bound 
state labeled by (12) of mass m(i2). If the mass of the bound state is 



3(12) 



(Re^ia^^ = 0, < ImO^^^^ < tt) 



1^(12) = v n^i + nT'2 + 2mim2 cosh 
the corresponding eigenvalue of the S- matrix Se{0) will have a pole at 6' = 6*12 such that 



Se{0) 



Ftp 



e - e^? 



for e^e 



(12) 
12 ) 



(2.10) 



giving rise to a residue Re- The eigenvalues are given by the diagonalization of the S- 
matrix 

^12(^) = E^e'^e(^)¥'l2 (2.11) 



where the projections onto the eigenspaces are given by the intertwiners (Clebsch-Gordan 
cofficients) ip\2 with 



12 - -^12 5 



p' 1 2 c- 
<^12<^e = 



Formula ( p.lOj ) may also be written as 

Res Si2(^) = 4'2)%2)V^Sf (2.12) 

""12 

where a matrix product with respect to the space of bound states V(i2) is assumed. 
Remark: In general an eigenvalue of S may have several poles corresponding to bound 
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states of different masses. On the other hand several eigenvalues may have poles at the 
same point, which means that there are several types of bound states (12) of the particles 
1 and 2 with the same mass. The space of the bound states V(i2) is then a direct sum of 
spaces belonging to these types of particles. 

The corresponding fields are related by a normal-product relation like 

^'(i2)(x) = A/'[^i^2](a;)</^;i% . 

The bound state S-matrix which describes the scattering of a bound state with another 
particle is given by ||68 



5'(12)3(6'(12)3) = \J R(12) Vl2^Slz (6*13) 5'23 (6'23 ) V^(12) / \J R(12) Q^_g^^gO-_^) (2-13) 

where the rapidity 6 (12) is fixed by pi + p2 = P(i2) • Here and below we use the phase 
convention that a/RI = %\/ —Re if Re < 0. 

In integrable quantum field theories there exist different types of bound state spectra 
which may be characterized by the absence or presence of solitons or kinks. Of course, in 
quantum field theory the bootstrap picture means that all particles are to be considered 
on the same footing. The names 'solitons', 'kinks' and 'breathers' are motivated by the 
classical non-linear equations associated with the quantum model. These equations may 
possess soliton or kink solutions, i.e localized non-singular solutions with a localized energy 
density. Special solutions consisting of a soliton and an antisoliton are called 'breathers' 
because of their oscillatory behaviour. In the quantum case we call a particle a soliton if 
it is a bound state of itself and another particle. Similarly (and more general), we call a 
particle a kink if it is a bound state of a particle with the same mass and another particle. 
The mass spectra of integrable quantum field theories characterized by the absence or 
presence of solitons are given as follows: 



i) There are particles labeled by a with mass 0, |6^ 



sm^ua 

ma = mi a = 1,2, . . . < 2/z/. 

sm 



This means that two particles of mass rria and form a bound state of mass me=a+b- 
The corresponding poles of the two particle S-matrix element and the rapidities in 
the bound state formula ( p.l3| ) are given by 



^ab = i^^io, + b) , 9a = 0e + i^ub , 6*6 = 6*^ - if z/o. 



The chiral S'f/(A^)-Gross-Neveu model |£2|, the Z{N) invariant Ising models 



or the SU{N)-a.Gine Toda field theories are examples for the above spectrum with 
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V = 2/N. In general the mass spectrum is more involved, for instance for affine 
Toda field theories (with real coupling constant) related to simply laced algebras 
the masses constitute the entries of the Perron- Frobenius eigenvector of the Cartan 
matrix |7I| and for theories related to non-simply laced algebras they do not 
even renormalise uniformly [|74 . 



ii) If there exist kinks (solitons) of mass M labeled by A then there are three types of 
bound states: 

a) Particles (breathers) labeled by a are kink-antikink bound states with 

m„ = 2Msin|i/a , a = 1, 2, . . . < 1/z/. (2.14) 

Here the corresponding poles of the kink-antikink S-matrix and the rapidities 
in the bound state formula (|2.13|) are given by 

d\j, = i7:{i-av) , eA = e, + \e\B , eB = e,-\e\j,. 

b) The kink B may be considered as a bound state of a particle a and a kink A 
such that the pole of the (a- A) -S-matrix and the rapidities in the bound state 
formula (|2.13| ) are 

OaA = y (1 + , = + y (1 - au) , Oa = 9b - iT^au. 

c) In addition as in i) two particles of mass rria and form a bound state of 
mass mc=a+b, however, here a < l/u. 

Examples for the latter case are the sine-Gordon (SU(2)-afiine Toda theory) alias the 
massive Thirring model with u = /^^/(Svr — /9^) = 71/(71 + 2g) and the 0(2A^)-Gross-Neveu 
model with u = 1/{N — 1). Also in this case the mass spectrum is in general more 
complicated, for example all affine Toda field theories with purely imaginary coupling fall 
into this category [|75| . 

f 12) 

The bound state formulae above may be depicted as follows: For 6*12 = 61 — 62 = O12 
with Im 612 > we introduce 

1(12) 21 11 

1(12) II 12 



1/^/^(12) {W)2? = r\ , <^(12)/\/%2) 

V o( 1^ ^ ^ V |^^2) . 
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Then we have the relations 




E 




where formally we have put Re/ Re = 1 even if i?e = in case that e does not correspond 
to a bound state. The sum in the last formula is over all eigenspaces Ve C Vi^ V2 of the 
S-matrix. If we would sum only over those e = (12) which correspond to bound states, we 
would get the projector onto the subspace of bound states in Vi ® V2. Moreover formula 
( p.l2| ) is depicted as 



Res 

(12) 



The bound state formula (|2.13|) may be depicted as 



(12) 






It implies relations of two-particle S-matrices called 'pentagon equations' (also referred 
to as bootstrap equations) like 



'S'(12)3 \JR(12) V\2^ = \[R(V2)V\2^ 81^82^ 





2.3 The Sine-Gordon model S-matrix 

The Sine-Gordon model alias the massive Thirring model is defined by the Lagrangians 

1 n 

respectively. 

The Fermi field correspond to the soliton and antisoliton and the bose field to the 
lowest 'breather' which is the lowest soliton antisoliton bound state. The precise relation 
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between the related coupling constants was found by Coleman []76[ within the framework 
of perturbation theory 



where the parameter u is introduced for later convenience. The two-particle S-matrix is 



/a 



S{9,u) 



\ 



b c 
c b 



S, 



bb 



\ 



(2.15) 



where the soliton-soliton amplitude a{9) and the soliton-antisoliton forward and backward 
amplitudes b{9) and c{9) 

a=X=X'^=X=X'^=X=X 



are given by 



b{9) 
a{9) 



sinh 9/1/ 



■a{9), c{9) 



sinh(i7r ~ 9)/v 

°° dt sinh 1(1 - z^)t 



exp 



sinh f, 

t sinh }^ut cosh |t m 



sinh ivr/z/ 
sinh(i7r — 9)/v 
9 



a{9) 



(2.16) 



These amplitudes fulfill 'crossing' 

aii'K -9) = b{9) , c{in -9) = c{9) 

and unitarity 

a{-9)a{9) = 1 , b{-9)b{9) + c{-9)c{9) = 1. 
The intertwiners iflf^ of section ^TI| are given by the non-vanishing components 



1, v>t, = yv2 

and the corresponding S-matrix eigenvalues are 

So = Sq = a , S± = b ± c . 



±1 



(2.17) 
(2.18) 

(2.19) 
(2.20) 



The amplitudes 5*0 = Sq have no poles corresponding to bound states. The amplitudes 
S±{9) have poles at 6' = in^l — ku) for even/odd k < l/u corresponding to the k-th 
breather as soliton-antisoliton bound states. 
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As examples of soliton-breather and breather-breather amphtudes those for the lowest 
breather are M 



, , sinh6' + i sin iyrfl + z/) dt cosh^z/t Q 

Ss^{B) = . • • n ^ =-exp / -2 -\-sinht- 

smh t7 — z sm ^vr^l + z/j Jo t cosh 2! «vr 

^ sinh 61 + 2 sin TTz/ r dt cos\i{\ — v)t 6 

Sbb{d) = —r^ — = - exp / — 2 ^ — smht— . 

smh tl — t sm nu Jo t cosh tir 



(2.21) 



The S-matrix element Sbb has been dicussed before in The pole of Ssb{0) at 9 = 

Z7r(l + z/)/2 belongs to the soliton as a soliton-breather bound state and the pole of Sbb{0) 
at 6 = mv to the second breather 62 as a breather-breather bound state. The intertwiners 
of section ^.1| are given by the non-vanishing components 



Vlb = ^Is = = 1- (2.22) 

The formulae involving higher breather may be found in e.g. 

^ , .L. dt cosh ^ut sinh ^ukt ^9 , , 

S,bM = (-l)^exp / -2 ' sinht- 2.23 

Jo t cosh Tjt smh T^z/r tir 



for k < I 



/■~ (it , cosh iz/tsinh iz/H cosh 1(1 - z//)t . , 9 

Sb,b,{9) = exp / -4 ? \ ■ ui I smht— 2.24 

Jo t cosh smh iz^t ztt 



and 



/■°° (it cosh iz/tsinh i(2A;z/ - l)A;t sinh i(l - z/)t ^9 , , 

Sb,bA9) = - exp / —2 2 sinht— . (2.25) 

JO t cosh 2^ smh 2^^^ zn 

3 Properties of generalized form factors 

We investigate the properties of generalized form factors, in particular for integrable 
quantum field theories in 1-1-1 dimensions. Some formulae, originally proposed in ||T^, are 
recalled and the physical arguments on how to derive them are provided in appendix |^. All 
arguments are solely based on the validity LSZ reduction formalism (see also [^) and 
the additional assumption of "maximal analyticity" which means, roughly speaking, that 
the S-matrix and the form factors are analytic functions everywhere except at those points 
where they posses singularities due to physical intermediate states. In other words the 
entire pole structure is of physical origin and in the following we investigate it employing 
the arguments of 0, 0, H (see also 0, [l|, Hf). 
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3.1 Form factors in momentum space and rapidity space 

For simplicity we first consider the case of bosonic cfiarged particles. Tfie extension to 
tlie general situation will be provided below. The corresponding Fock space is spanned by 
the in- or out-states of particles and anti-particles given by ( |2.1| ) and (|2.2D. In addition 



to the notation of vectors and matrices of ( p.6| ) we denote co-vectors by 

Vl...n e (3.1) 

with components Va = Vai„,a„- 

Let now 0{x) be a local scalar operator, the generalized form factors are defined as 
the co-vector valued functions given by 

(0 I 0{x) . . .,Pn)Z...c.„ = e-(P^+-+^") F^{{s,, + ze)(i<.<,<„)) (3.2) 

where a = {ai, . . . , a„} and where = (Pi+Pj)'^ is one of the Mandelstam variables, as in 
the previous section. There may also be anti-particles in the state. As is well known these 
functions are boundary values of analytic functions as indicated by the e-prescription. We 
assume that the domain of analyticity is much larger than could be proven by means of 
general principles. Similar as for the scattering matrix we assume in addition at this 
point "maximal analyticity" meaning that there should be no redundant poles, but all 
singularities should be of physical origin as particle states etc. Since the x-dependence of 
the form factors is trivial, in the sense that we may always carry out a translation as in 
eq. (|3.2| ), we consider in the following the operator always at the origin, i.e. O = O{0). 

Under the assumption that F is an analytic function, an interchange in eq. ( |3.2| ) of 
the in and out states leads to the replacement of s + ie by s — ie. This means in particular 
that 

( I O I pi, . . . ,p„ = F^{{S,, - Z6) (!<,<,<„)) . (3.3) 

The crossing property for the connected part of the matrix element yields 

a,..Z {pi,...,Pm\0 I p^+i, . . . , p„ = ^^i^^J + ' ^ki + ^e) (3.4) 

where 1 < i < j < m, l<r<m<s<n, m < k < I < n and t^s = {Pr ~ PsY another 
Mandelstam variable. See appendix ^ for a proper derivation of this claim. 

The most basic properties of the form factors are usually refered to as Watson's equa- 
tions |]7^, which have been already known in the fifties. It is instructive at this point to 
discuss them first for the case n = 2. Using the completeness of the out-states we have 

FZ.2 (^12 + 26) = ( I O I pi , P2 a. =Y.{^\0\0Ut){ OUt \ p,, p, „^ . (3.5) 

out 
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For 4m^ < S12 < 'lowest inelastic threshold' only the two particle S-matrix contributes 

+ ^e) = ^^5.^,(^12 - Stltlis,,) (3.6) 

and analogously starting with ai (pi | | )a2 in ( P-5|) we obtain 

FZc.,{ti2 - ^e) = i^S„,(ti2 + le) , (3.7) 

where the fact has been used that the one-particle S-matrix is always trivial. In integrable 
theories there are no inelastic transitions, therefore eq. (|3.6|) holds for all s > 4m^. The 
generalized Watson's equations for 1 < m < n read as (see |]13|) 

Faisij+ie, trs-ie, Ski+ie) = 5'^r;.„?^'(sij) F^,{sij-ie, Us+ie, Ski-ie) 5^1+ (s/^). (3.8) 

For a diagonal S-matrix these equations have been discussed before in []75| . 

The generalized form factors also contain singularities |T3|, |T^, |^ which are determined 
by the one-particle states in all sub-channels (a^, . . . , a^) C (ai, . . . , a„) (see figure 0). 
Poles occur if the square of the total momentum in the sub-channel equals the one-particle 




Figure 2: A singular contribution to the n-particle form factor diagram corresponding to 
a sub-channel. The dashed lines belong to off-shell lines. 

mass squared. In particular there are poles, if for instance particle 1 is the anti-particle of 
particle 2 and particle 1 is crossed to the out-state together with p2 ^ Pi, which means 
(P2 + Ps — Pi)'^ f^l- Alternatively, if particle 3 is a bound state of particle 1 and 2, in 
which case (p2 + Pi)^ "^3- The residues of the form factors at these poles are related 
to form factors with fewer legs, as indicated in figure ^ We will discuss these facts later 
in detail. 

Similarly as for the S-matrix we may also write the form factors (|3.4| ) as co-vector 
valued analytic functions of the rapidity differences % = 6i — 6j 

F^{sij + ie,trs-ie,Ski + ie) = F^{\eij\,i7r - \9rs\,\^ki\) 
F^{sij - ie,trs + ie,Ski - ie) = F^{-\eij\,iTr + \9rs\, -\Oki\) 

The domains of analyticity and the physical regimes in the complex planes of the Man- 
delstam variables and the rapidity difference variables are depicted in figure ^ However, 
now the branch cut between between region II and IV is absent (c.f. eqs. (|3?^ and (|3.8|) ). 
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3.2 The auxiliary form factor function 



' a \ 



Furthermore, it is convenient to introduce a new co- vector valued auxiliary function 
which is considered as an analytic function of the individual rapidities of the particles, 
instead of analytic functions of all rapidity differences (see also []T7[)- It coincides with 
the generalized form factor for a particular order of the rapidities 



/f(^i,...,^„) = F„^(|%|) = (O|0|pi,...,j9„)r , for ^i>...>^„. 



(3.9) 



For all other arrangements of the rapidities the functions f^iO) are given by analytic 
continuation. The domains of analyticity, the physical regimes and the transitions to the 
crossed regions in the complex planes of 6i and 6j for 6i > 9j are depicted in figure ^. 



+ 



II 



II 



+ 



III 



II 



(a) (b) 

Figure 3: The physical regimes in the complex planes of the rapidity variables (a) 6i and 
(6) 6j for 9i > 9k > 9j {k = 1, . . . ,n) . Again the crossing transitions (see appendix B) 
are indicated by the arrows. 

Now we formulate the main properties of generalized form factors in terms of the aux- 
iliary functions /f „ under the assumptions of "maximal analyticity". 

Properties: The co- vector valued auxiliary function /fj' „(^) is meromorphic in all vari- 
ables 9i, . . . ,9n and 

(i) fulfills the symmetry property under the permutation of both, the variables 9i,9j 
and the spaces i,j at the same time 



f...ij..{- ■ ■ ^j, ■ ■ ■) = f...ji..{- ■ ■ 1 ^j, ^i, ■ ■ •) Sij{9i - 9j 
for all possible arrangements of the 6''s, 



(3.10) 



(ii) fulfills the periodicity property under the cyclic permutation of the rapidity variables 
and spaces 

A^..„(^i, 92,...,9n,) = fZ.ni{02, ...,9^,9,- 2m), (3.11) 
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(iii) and has poles determined by one-particle states in each sub-channel (see figure 0). 
In particular the function fa{§) has a pole at 612 = in such that 

Res fZ.niOl, . . . , = 22 C12 fZ.ni03, ...,9n)(l-S2n... S23) (3.12) 

where C12 is the charge conjugation matrix with matrix elements Caa' = ^aa'- 

(iv) If there are also bound states in the model the function fa{0 has additional poles. 
If for instance the particles 1 and 2 form a bound state (12), there is a pole at 



6*12 = Orz'^ such that 



f^eS /i^...„(6'i, 6^2, ... , 9n) V3(i2) = f {12)... ni.^ {12), ■ ■ ■ , 0^) ^j2iR(i2) 



o 



(3.13) 



where the quantities v^(i2); -^(12) and the values of 6*1, 62, 0{i2) and 612 were dis 



cussed in section 2.2 



The property (i) - (iv) may be depicted as 
(i) 



/ 



iii) 



[III 



{iv) 



— Res 

2i ei2=j7r 



Res 





) - 


















n 



/ 



o 




f 



o 



Both properties (ii) and (iii) are consequences of the general crossing formulae 

l{pi\0\p2,...,Pn)t.n 

n 

= E I Pj )j f2...j...n ^'^i ■ ■ ■ + /S...n(^i + «vr„, . . . , On) (^3 J4) 

n 

= E i{Pi\p,),f^- ^S,^---S,,+, + fl,^,{...,e^,9^-in.)C'\ 

where we introduced the notation ^, meaning that the space j and the corresponding 
variable 9j are missing. In terms of the components, i{pi \ pj )j means ^q^qj 47r5(6'i — 9j) 
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and S^^ means Sa^ai- These are equations for distributions where on the right hand side 
the second terms are understood as boundary values of analytic functions with 7r_ = tt — e. 
The crossing formulae may be depicted as 



/ 




f 



o 




f 



o 



where we have again used the graphical rule (|2.9|) , which states that a line changing the 
"time direction" also interchanges particles and anti-particles and changes 6 6 ± in. 
Taking the analytic part of the crossing relation one obtains property (ii) and considering 
in addition the part with point like support one gets property (iii). The proofs of the 
properties (i)-(iv) and equation ( |3.14|) are provided in appendix 0. 



(v) Naturally, since we are dealing with relativistic quantum field theories we finally 
have 

(3.15) 
e'^^'^O where 



A^..jei + n,...,^„ + M) = e^^"/L..„(^i, 



±SU rO 



if the local operator transforms under Lorentz transformations as O 
s is the "spin" of O. 



3.3 The general bosonic and fermionic case 

For the general case where the states involve also fermions and where 0{x) is a local 
bosonic or fermionic operator with arbitrary spin we write the matrix elements of O{0) 

as 



( 1 o I Pi, . . . , P. )Z...^„ = Y.^iPi) ■ ■ ■ rS...,, ■ ■ ■ u{p,) ■ ■ ■ p^: g£) 



Sij + ie) (3.16) 



where the F are matrices in spinor space. For the invariant form factor functions G^, 
the Watson's equations look quite analogously to those in the bosonic case. However, 
sometimes it is more convenient to consider the full matrix elements and then we must 



take into account sign factors due to the fermions. Analogously to eq. (|3.9|) we introduce 
the CO- vector valued auxiliary function which determines the form factors for a specific 
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order of the rapidities. For the general case the three main properties of the co-vector 
valued function may be written as: 

(i) fL..ij...n(^l) ■ ■ ■ ■ ■ ■ ,(^n) = fl...ji...n{(^l-, ■ ■ ■ , Oj , 6i, . . . , 6.^) Sij 

(iii) ^ ^Cufl.M,---,0n){l-S,^...S,,) 

(3.17) 

The bound state formula (iv) is in general true for the invariant part of the form factors. 
For the case of fermions, spinors have to be taken into account (see the examples below). 
In the formulae ( |3.17| ) the statistics of the operator O is taken into account by aoi = — 1 if 
both O and particle 1 are fermionic and aoi = 1 otherwise. The statistics of the particles 
is taken into account by S which means that S12 = —Su if both particles are fermions 
and S12 = S12 otherwise. Again, both properties (ii) and (iii) are consequences of the 
crossing formulae, which, for the general case of bosons or fermions, reads 

l{pi\0\p2,...,Pn)t.n 

= <Toi I I ), /?..-,..„ ^2, ■ ■ ■ 4-1. + G''f?2..A0i + • • • , ^n)| (3.18) 

= E -. „^,„---%+i + /2^..„i(...,^„,^i-^vr_)Cii. 

j=2 



replacing eq. ( 3.14 ). The proof of these relations are also given in appendix 



The appearance of S is natural in the context of factorizing S-matrices. See for example 
the general Yang-Baxter relation ( |2.5| ) which is essential if transitions as fermion + anti- 
fermion — »• boson -|- anti-boson are possible. 

4 Solution for the sine-Gordon alias massive Thirring 
model Model 

We will now provide a constructive and systematic way of how to solve the properties 
i)-v) for the co-vector valued function / once the scattering matrix is given. To capture 
the vectorial structure of the form factors we will employ the techniques of the "off-shell" 
Bethe ansatz which we now explain briefly. 

4.1 The general formula 

As usual in the context of algebraic Bethe ansatz we define the monodromy matrix 



Tl...nfi{d_, 6q) = Siq{6i — 6q) 320(62 — do) ■ ■ ■ Snoi^n — ^0 











1 


2 


n 






(4.1) 
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as a matrix acting in the tensor product of the "quantum space" Vi...„ = Vi ® ■ ■ ■ ® 
and the "auxihary space" Vq (all Vi = C'^ = soliton-antisoliton space). The Yang-Baxter 
algebra relations yield 

Tl...n,aid, Go) Ti,,,n,bid, Oh) Sabi^a — ^b) = Sahi^a — Oh) Ti,„n,b{0, Ob) Ti^..n,a{0_, Oa) (4.2) 

which in turn implies the basic algebraic properties of the sub-matrices A, B, C, D with 
respect to the auxihary space defined by 



Al...n%0) fil...n(£( 
Cl...n{0.0) D,,„n{0,i 



(4.3) 



A Bethe ansatz co-vector in Vi ^ is defined by 

'4^l...n{0,Ui, . . . ,Um) = ^l...nCl...n{0,Ui) ■ ■ -Ci^niOyUm) 



Or, 



Or, 



Ul 



(4.4) 



Ur 



where f^i...n is the "pseudo-vacuum" co-vector consisting only of particles of highest 
weight. When the monodromy matrix involves only the scattering matrix of soliton anti- 
solitons it is given as 

J^l...n =T ® ■ ■ ■ ® T (4.5) 

consisting only of solitons and fulfilling 

^l...nBi,„r,{0,u) = 

n 

^l...nAi,„niO,u) = Y[a{Oi- U)ni,„n 



i=l 
n 



(4.6) 



ni...nD,,„^{0,u) = l[b{0,-u)n,. 



1=1 



Here the eigenvalues of the matrices A and D, i.e. a and b are related to the amplitudes 
of the scattering matrix (refer ( p. 161 ) via a = —a and b = —b. 



In the conventional Bethe ansatz [£0|, one is usually concerned with the computation 
of the eigenvalues of the transfer matrix 

ri...niO, u) = Ai...„(^, u) + £>!...„(£ u) 

on a Bethe wave vector. Applying the transfer matrix to the co- vector 
in general an equation of the form 

'ipl...n{0, Ul,. .., Urn) ^1. ..„(£, u) = A{u\ui, . . . \0) i)l...n{0, Ml, ... , Um) 



(4.7) 
one obtains 



(4.^ 



- ^ Aj(ni, . . .,Um\&) i'l^AQ-Wu . . . , Mj-l, M, Uj+i 



, . . . , Um 
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where A{u\ui, ■ ■ -10) and Aj{ui, . . . ,Um\0) are some complex valued functions. The co- 
vectors 4'i„,n{i.\ui, . . . , Uj-i, u, Uj+i, . . . , Um) are not proportional to the Bethe ansatz 
vectors ipi...n{0,Ui, . . . ,Um)- Hence in general, that is for an arbitrary set of the spec- 
tral parameter, the Bethe ansatz vector is not an eigenvector of the transfer matrix. 
To achieve this one usually imposes the validity of the Bethe ansatz equations, i.e. 
Aj{ui, . . . ,Um\0) = (j = 1, ...,m) such that the so-called "unwanted terms" vanish 
and one obtains a genuine eigenvalue equation for the transfer matrix with eigenvalue 
A{u\ui, . . . \6_). In analogy to the one particle situation one may refer to such Bethe vec- 
tors as being "on-shell" in contrast to the generic situation ( [4.8| ) which is referred to as 



"off-shell" 0. In order to construct solutions to the properties (i)-(v) we shall employ 
the Bethe vector ( |4.4| ) in its "off-shell" version. 

Let us now consider the auxiliary form factor function given by 

f^{9,,...,9n) = {0\O\pi,...,pny^ , for 9,>...>9r.. (4.9) 
where the indices a refer to solitons and antisolitons. 

Theorem 4.1 The co-vector valued function /f „(^) fulfills the conditions (i), (ii) and 
(iii) of section ^ (see eqs. (\3.1(\\3.1^ )) if it is represented by the following generalized 



Bethe ansatz / |5q/ 

f?..M)=N^ I du,---l dUrng{9,u)Q,...nC,...n{9,Ul)---Cl...n{i,^^rn) (4.10) 



'Cg_ JCo_ 

with a normalization constant and the scalar function 



n m / „ „ \ 

n m,)i[i[m-u,) n r(«.-«,)e^^'^"^-^'^J (4.ii) 



l<i<j<n i=l j=l 1<'<J<"^ 



where s = s/q and s is the "spin" (c.f. eg. ( \3.13i )) and q = n — 2m is the charge of the 
operator O. The number s is assumed to fulfil exp(2'7ris) = (—1)"- The function F{9) 



(see ([j.lJl )) is a soliton-soliton form factor fulfilling Watson^s equations 

F{9) = -F{-9)a{9)=F{2m-9) (4.12) 



with the soliton-soliton scattering amplitude a(9) (see \2.1d[ ). The scalar functions (j){u) 
and t{u) are defined as 

4>{u) = ^, X r}/ ^ ; '^{'U) = ,/ N ''^z r • (4-13) 

^ ^ F{u)F{u + t7T) ' ^ ^ (j){u)(f){-u) ^ ' 

The integration contour Ce consists of several pieces (see figure^): 
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a) A line from — oo to oo avoiding all poles such that Q6i — tt — e < '^uj < '^di — tt. 

b) Clock wise oriented circles around all poles (of the (j){9i — uj) ) at uj = 6i. 
In addition we assume that the number of particles involved, i.e. n, to be odd. 



o ^1 +m(2z/- 1) 

o 9n + in{2u - 1) o 02 + ^vr(2z/ - 1) 

• 6i + in 

• 9n + iiT •62 + 111 

o di + «vr(z/ - 1) 

o^„ + ot(z/-1) o ^2 + ^7r(z/ - 1) 



o 6^2 - ?7r f~ 



» On — 2TTi • 62 — '^TCi 






— m 




- i-RV 


•01 


- 2m 



Figure 4: The integration contour Cg (for the repulsive case v > 1). The bullets belong to 
poles of the integrand resulting from a{6i — Uj) (piOi — Uj) and the small open circles belong 
to poles originating from h{9i — Uj) and c{6i — Uj). 

This theorem is proven in appendix |B[ 

Remarks: 

• The minus sign in eq. (|4.12| ) is due to the fermionic statistics of the sohtons. 



A solution of the Watson's equations ( [4.12|) is 



F{9) = -i sinh fZ'^ie) (4.14) 

where the 'minimal' soliton-soliton form factor function is given as 

r^^ie) = exD r sinh 1(1 -z/)t l-cosht(l-g/(27r)) 
^' Jo t sinh iz/t cosh |t 2 sinh t 

The corresponding functions 0(m) and r(-u) are (see appendix |C|) 

1 dt sinh |(1 — i')t ( coshi(:(i _ (in)) — 1 



[u] = const. — - — exp / 
sinh u Jo 



t sinh I ut sinh t 

t{u) = const, sinh M sinh m/z/ 
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Using Watson's equations ( [4.12| ) for F{u), crossing ( |2.17| ) and unitarity ( |2.18| ) for the 



sine-Gordon amplitudes one derives the following identities for the scalar functions 
(f){u) and t{u) 

,, , ,,. . 1 ,/ • \ aiu — 2m) /. 

(piu) = (pnn — U) = —-r-TT (pyu — ITT) = -— (plu — ZTTt) , (4-15) 

b{u) h[u) 

/ N Hu] a(27Ti — u) , , , , 

r(u) = r(-u) = ^ -j frn-27rz 4.16 

^ ^ ^ ^ a{u) b{2TTi -u) ^ ' ^ ' 

where h{u) is the soliton-antisoliton scattering amplitude related to a{u) by crossing 
h{u) = a{m — u). 

The number of C-operators m depends on the charge q = n ~ 2m of the operator 
O, e.g. m = (n — l)/2 for the soliton field il){x) with charge (? = 1. 



The integrals in eq. ( |4.11| ) converge if |(1 + l/'^)^' =f2s + 2/z/+1 >0. 



Note that other sine-Gordon form factors can be calculated from the general formula 
( [4.10| ) using the bound state formula ( |3.13| ). 



We shall now apply the general formula (^.1CI|) to an explicit example and exploit the fact 



that the properties (i)-(iv) relate several different form factors to each other. This will 
permit us to carry out various consistency checks. 

4.2 The two particle form factors 

We repeat some well known results (see for example p!3| , |69[]). According to equation 
( p. 17] ), the auxiliary function for the two particle form factor /^(^i, O2) has to satisfy 



a' 13' 

These matrix equations may be solved easily by diagonalization of the S-matrix. If there 
are only bosons involved we have to solve the scalar "Watson's equations" 

ff{0) = f?{-0) Se{0) = ff{2m - 6) (4.17) 



where Se{0) are the eigenvalues of the S-matrix given by eq. ( |2.11| ). In |J13[ it was shown 
that the general solution of these equations is of the form 



f^{0) = N^K,{0)fr''{0) (4.18) 
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where Nf is a normalization factor, f^^"'{9) is the minimal solution of Watson's equations 
without any poles or zeroes in the physical strip < Im ^ < tt and Kg{6) is an even periodic 
function with period 27ri. If the S-matrix eigenvalue is given by 

roo 

Se{e)=exp dtf{t)smhte/in (4.19) 
Jo 

the minimal solution of Watson's equations is given as 

1 - cosht(l - e/i-K) 



exp / dt f (t) 
h 



2 sinh if: 



(4.20) 



If there are also fermions involved Watson's equations (|4.17| ) hold for the invariant form 
factors (c.f. (|3.16|) ). For the full matrix elements the representation (|4.18|) holds with 
additional factors exp(±^^j/2) on the right hand side for all fermions. 

The poles of f^{0) in the physical strip are determined by the one-particle states in 
the channel corresponding to the S-matrix eigenvalue. In [|13| the minimality assumption 
was made, meaning that there are only these poles and no zeroes in < ^3^^ < vr. This 
implies that 



Ke{e) = n 



1 



k=i sinh I ( 



9k) sinh ; 



0, 0<'^e.<7r) 



For several examples this assumption was checked against perturbation theory. 



Examples 

We present two-particle form factors for several local operators and several particle states 



of the sine-Gordon quantum field theory. Some of them were already calculated in |13 



see 



also fl^)- Up to normalizations the problem is solved by eqs. ( |4.18| - |420|) since the sine- 



Gordon S-matrix (|2.15|) is diagonal except of the soliton-anti-sector where the eigenvalues 
are given by eqs. (|2.16| ) and (|2.20|) . 



4.2.1 The two-breather form factor 

The simplest sine-Gordon form factor is that for a scalar operator (j)'^{x) = M(j)'^{x) con- 
necting the two-particle lowest breather state to the vacuum 

(^12) = ( 1 0^ I = Ni K^eu) /r"(^i2) . 



According to ( |4.19| ) and ( [4.20| ) the minimal form factor function combined with ( |2.21| ) 
reads 



fl^'^iiiTx) = -zsinhi^ exp 



°o dt cosh(i - iy)t 1 - cosht(l 



x] 



t 



cosh 



2 sinh t 
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The "minimality assumption" implies that the 'pole function' 

^''^^^ " sinh 1(0 - ivrz/) sinh + ZTTZ/) 

only possesses the pole corresponding to the second breather 62 as a bound state of 
two lowest breathers b. The normalization constant can be calculated by means of the 
asymptotic behavior |T^. Weinberg's power counting implies that in the limit of infinite 
momentum transfer the form factor tends to its free value 

{0\<P'\puP2)T,^{0\:^':\puP2)lr = 2Z^ as (p, + p,f ^ -oc . 
Here Z'^ is the wave-function renormalization constant of the fundamental sine-Gordon 



field which has been calculated in \13 



Z^ = {l + v)-^E{v). (4.21) 



sm 



1 j-KX + 

E{x) = exp ( / - — dt ] (4.22) 



We introduce the function 

71 Jo sin if: 

= exp I - i (^2 L (e*'^^') + t L [e'^'^ + l) - x In (e''^" + l) + | 
00 

where L{x) = ^ + | lnxln(l— x) denotes the Rogers dilogarithm ^0|- Notice, that this 

n=l 

wave function renormalization constant satisfies (compare figure ^ < Z*^ < 1, which is 
a general consequence of positivity |^ (see also e.g. [0] p. 204). For the free boson case, 
that is /3 = 0, u = 0, we have Z"^ = 1. For the free soliton case, i.e. (7 = 0, z/ = 1, where 
the breather decays into soliton-antisoliton pairs, we have Z'^ = 0. Using the asymptotic 



formula the normalization has been calculated in 13 



iV,f = -2(l + z/)^z/cot|z/. 

In addition we may now employ this result and compute a further renormalization 
constant by means of the bound state formula (p.l3|) . With ipl^^ = 1 (see eq. ( p.22|) ) we 
calculate the wave function renormalization constant Z'^' via 

Res (6) {2z Res SM) = f£ = ( | 0^ | Pi+P2)Z = ^ 

where the fusing angle is 9q = mv. The wave function renormalization constant turns out 
to be 
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Figure 5: The wave function renormalization constants Z"^ and Z*^^ as a function of the 
couplings =^ = 

Aerain we have < Z"^ < 1 and now Z*^ vanishes at u = and u = 1/2, which is to 
be expected since at these values (compare for instance with the mass formula ( p.l4| )) 



the second breather decays into the two lowest breathers or a soliton-antisoliton pair, 
respectively. 

4.2.2 The breather-soliton form factor 

We now choose 0{x) to be the fermi- field ip{x) of the massive Thirring model which 
annihilates the soliton. We assume the breather-soliton form factor related to this field 

ft{e,,e2) = {o\ij\PuP2)Ts 

to acquire the following form 

ftie^, e^) = Nt (l + Nf f coth leu)u{p2) KUOu) /r"(^i2) , (4.24) 

that is consisting out of a scalar and a pseudoscalar coupling part. For our conventions 
concerning spinors and the 7-matrices see section 5. Upon employing ( [4.19| ) and ( [4.20|) , 
the minimal form factor function reads together with equation ( |2.21| ) 

rminr \ ■ n /"^ rft COsh ^Z/t 1 - COsh t (1 - x) 

U. (™.)=smf.exp/^ T'^Sh|r 2^t ' 

Extracting explicitly the expected pole structure the "pole function" reads 

^''^^^ ^ sinhi(^-^7ri±^) sinhi(e + mi±^) ' 



27 



Once more we may use the bound state formula ( p.l3| ) in order to compute the normal- 
ization constants, with (p^f = 1 and 9o = i'^{v + 1) (see eq. (|2.22|)) to obtain 



Res /r.(^i,^2) (2z Res 8^612, 

«12=6'o ^ 012=00 



2 

-1/2 



/J = (0|^|pi+P2)r = «(Pl+P2) 

This determines the normalization constant to be 



1 TT 

COS -^1/ 



hs 



and the ratio of pseudo-scalar and scalar coupling to be 



Nf 



tan |z/ tan |(1 + z/). 



Note that formula (WM) may alternatively also be written as 



Nt cos^^fl-z/) 



^ sinf(l + z/ 



liPl +P2) -rn 

_ 1 
e 2 



1 + 7' COthi^l2)«(j92)/r"(^^12) 



J ITT 1/7^ 



+ 



«(P2)/r"(^ 



12; 



cos 



2 f z/ I sinh i(^i2 + ^0) sinh 1(^12 - 9, 



sinh \9i2 



4.2.3 Soliton-antisoliton form factors 

Having an in-state, which involves a soliton and an antisoliton, we have several options 
for the operator 0{x) such that the two-partice form factor is non vanishing. 



a) Let 0{x) 



N'lp'-f^tplx) the electromagnetic current [|13 



The function f-{9) fulfills Watsons equations with the negative C-and P-parity 

- Ihir^li/X ^) (seeeqs- mMM)- Taking the 



S-matrix eigenvalue S^{9) 



singularity structure into account we obtain [|T3|, with the help of ( [4.19|) , ( [4.20|) and 



and 



J S. 



niTx] 



f-io) 



exp 



cosh j{in — 9) 
cosh ^{in — 9) / u 



dt sinh i(l - z/)t 1 - cosht(l 



x) 



t sinh iz/t cosh 



2 sinh t 
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b) Let 0{x) = (f){x) the fundamental sine-Gordon field which correspond to the lowest 
breather UTT 



sinh 6*12 



f- 



712 J 



The function f-{9) is the same as in a). Since Coleman's correspondence |j76| relates 
the field (p and the current j'^ by 



27r 



(4.25) 



the normalization constant turns out to be [|T^ 

27ri 



Ni 



We may now carry out a consistency check and compute once more the wave function 
renormalization constant now starting, however, from a different form factor. For 
this purpose we use once again the bound state formula (|3.13|) with ipl'^ = —^py = 
l/\/2 (see eq. (|2.19|) ) and calculate 

Res ff,{e) V2(2i Res S-{e)^ '^''^ 



Oo 



f^={0\<l>\p)r = VZ^ 

where the fusing angle is 6*0 = ^{1 — v). This computation leads to the value for 
wave function renormalization constant of the previous subsection ( |4.21| ) which has 
been obtained in []13[ by slightly different arguments. 



c) Let 0{x) 



{ 1 0' I PUP2 )s = Nf; mnoi) u{0i2) ■ 



The function f^{6) fulfills Watsons equations with the positive C-and P-parity S- 
matrix eigenvalue S+{e) = -'^^i^^zj^ ci{e) (see eqs. (^JWM )- With (|4l9|) 
and ( [4.201 ) we obatin ||l3l together with the explicit expression for the integral rep- 
resentation of this amplitude of the scattering matrix (|2.16|) 



sinh |(i7r — 9) 
sinh |(i7r — O)/^ 



/—(e) . 



We still have to fix the normalization constant Nf^ , which may be achieved by 
employing the bound state formula ( |3.13|) . Taking ipl^^ = Lfl^ = l/\/2 (see eq. ( p.l9| )) 
we calculate the wave function renormalization constant to be 

1/2 



Res ft; (9) V2 {2t Res 8^(9)) " = /,f = ( | 0^ | p )S = 
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where Oq = i7r(l — 2z/). The wave function renormahzation constant Z"^^ was cal- 
culated in the previous subsection ( [4.23 ), such that we obtain the normalization 



constant 

8M sin^ f z/ ■ 

4.3 Three particle form factors 

We shall now analyse the general expression proposed in theorem 4.1 for an explicit 
example. First we recall the three breather form factor which was already calculated in 
p!3| and apply a consistency checks using (iii). Furthermore we calculate the three soliton 
form factor using the general formula ( [4.1(j| ) and apply some consistency checks using (iii) 
and (iv). 

4.3.1 The three breather form factor 

We choose the operator 0{x) to be the fundamental sine-Gordon field (pix) which corre- 
sponds to the lowest breather and consider the form factor 

fLiOue2,es) = {o\<p\puP2,P3)Z- 

The minimality assumption suggests the proposal 

fLiOi, ^2, ^3) = ivi K,uoi, ^2, ^3) fnou) /r"(^i3) /r'^(^23) 

with the "pole function" 

= coshi... coshV.3 cosh I.., ■ 

The "two-breather pole function" Khi,{6) and the minimal form factor /^™(^) were already 
provided above. We use property (iii), i.e. the recursion relation ( |3.12D , and calculate 



Res /l(^i2, ^13, ^23) = 2tVZ<t> (1 - Shb{e23)) 



which determines the normalization constant ||T^ 



Nl^^j^ = + uf cot f z/ cos^ f z/ (Z^ 



-3/2 



30 



4.3.2 The three (anti)-soHton form factor 

We now choose 0{x) to be the fermi-field ip^lx) of the massive Thirring model which 
annihilates the soliton. We consider the form factor 



/?23 (^1,^2,^3) = (0 I V^^|pi,P2,P3) 



123 



Here ± refers to the first or second component of the spinor, respectively. Nonvanishing 
matrix elements contain two solitons and one antisoliton. Taking in the general formula 
( [4.10| ) for n = 3 and m = 1 we obtain 

3 

ft2-M) = Nt n Fie,,) [ duY[m-u)e^^^''~^''^'^ni23Ci23{e,u). (4.26) 

i<«<i<3 1=1 

Here the function F{6), which fulfills Watsons equations 

F{e) = -F{-e) a{e) = F{2'ni - e) , 

is closely related to the minimal form factor which was computed above 

F(^) = -^sinhi^^/r"W. 

The scalar function reads 



F{u)F{u + m) 

We now use property (iii), i.e. the recursion relation ( |3.12|) and calculate 



Res ft^M = 2^C,2 ft (1 - ^23(^^23)) 



which determines the normalization constant 



Nf = ±t^{frr{0)y. (4.27) 

Note that this follows also from the general recursion relation ( p. 7]) . The form factor is 
now fixed with all its constants. However, we also expect the bound state formula ( |3.13| ) 
to hold and we may employ it now as a consistency check. We calculate with (p^ given 
by eq. ( |2.19|) and the fusing angle given hj Oq = i7i{l — u) 

e^^ta^ a /j ^ ,.12 /o^- c c (a\\'^'^ - fi'^^ 



012=00 ^ 0=00 ' 

The result of this computation coincides with the form factor proposed in ( [4.24[ ). Having 
convinced ourselves of the mutual consistency of several solutions we shall now carry out 
an additional check and compare the results with conventional perturbation theory. 
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5 Perturbation theory: Massive Thirring model 

In order to check the three particle form factor of the fundamental fermi field correspond- 
ing to the soliton in perturbation theory we calculate the four point vertex function. We 
start with the Lagrangian 

The fermi field i\}[x) annihilates a soliton and creates an antisoliton with the following 
normalisation 

( I ^(x) I p = e-*P^ u{v) , ( I ^(x) \p)^ = 6^, e'^^ v{p). (5.28) 
We use the following conventions for the 7-matrices 

and for the spinors 

uip) = ( '^"//f ) , v{p) = v^z ( ^"J/' ) with = M ( ^^^^l^ . (5.30) 
We also employ the formulae 

{7^ Y} = '^g'" , [Y, 7l = ^e^'^Y , {e^'' = -e''^ , = 1) 
7'^7'^7^ = , e'^^e^'^ = g'^'^g''^ - g^'^g^'^ , = e^^^p . 

The Lagrangian implies the Feynman rules of figure |] 

= -ig-f>' ® 7^ , J! - 




-fk- M 

Figure 6: The Feynman rules for the massive Thiring model. 
The three particle matrix element of the fermi field up to order g turns out to be 

0\m\Pl,P2,P3)Tss = 

Pi P2 P3 

= -ig —, \ \ ^ 77 {l^u{p2) v{pi)^f,u{p3) - YuiPs) v{pi)^^u{p2)) + 0{g^) 

l[Pi +P2+p-i)- M ^ / 

- -ig sinh ^9 ^(P^) cosh \9i2 + u{p3) cosh 16*13 ^ ^,2-. 
^ coshi^i2 coshi6'i3 coshi6'23 
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Note that for the sohton-sohton scattering amphtude this imphes 



a{9) = l-ig tanhy + 0{g^ 



in agreement with eq. ( p.l6D 



To calculate the exact form factor up to this order we start from the general formula 

(H) 

/£(^l,^2,^3)= / dul{9,u) 



with the integrand 



l<«<i<3 i=l 



u 



ss 



U) 



Using the residue theorem the integral may be written as 

/ dul(6,u) = 27ri I Res -if Res + Res + Res- Res ) ) Ji 

7Ce \^6li-j7r ^^1 62 e-i 6»i+i7r(i.-l) / y 

+ \ I du(i{e,u) + i{e,u + i7c)) 

JCo ^ 

where the integration contour Co is a line from — oo to oo avoiding all poles such that 
+ Tc{iy -2) <^u< (for u > 1, u = 1/(1 + 2g/Ti) ^ 1). The integral on the right 
hand side is of higher order in g and the residues give 

f^^ff) f) /3 ^ _ /V'^^^l!!^ Ifl n^(p2) cosh igi2 + n^(p3) cosh |gi3 n^ 

VM cosh Tj6'i2 cosh 7^6^13 cosh 7^623 



which is consistent with the result of the Feynman graph calculation because of equation 
( [4.27| ) and /™"(0) = 1 + 0(5'). Hence we obtain mutual consistency between the solutions 



of the form factor equations and conventional perturbation theory. 



6 Conclusions 

We have outlined in detail the so-called "form factor program" . Using only the "maximal 
analyticity assumption" and the validity of the LSZ formalism we have derived general 
properties of form factors. The properties are expressed in terms of the equations (i)-(v). 
We provide a solution for these equations in a closed form, which captures the vectorial 
structure by means of the "off-shell" Bethe ansatz and the singularity structure in term of 
certain contour integrals. The validity of this solution has been checked by constructing 
various explicit two and three particle form factors. We have compared our solution for 
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the three particle form factor of the fundamental fermi field with the expressions obtained 
from perturbation theory in the massive Thirring model. We find complete aggrement 
between these two approaches and we take this as a further indication for the validity of 
the "form factor program" formalism. 

The vectorial nature of the form factors we present, is encapsulated in the "off-shell" 
Bethe ansatz and the singularities are encoded in certain contour integrals. We assume 
that this structure is of a universal nature and will allow to construct further solutions of 
other integrable theories. It will be highly interesting to work out such solutions explicitly. 
This task and the detailed study of the correlation functions obtained from these solutions 
is left to future investigations. 

We have applied the general formula ( [4.10| ) to an explicit example and exploit the 
fact that the properties (i)-(iv) relate several different form factors to each other. This 
permits us to carry out various consistency checks. We have for instance the following 
relations 

Jbbb ^ Jb ~ h J ssb ^ Jbb 



rip y-") 
J b2S 


rip (i^) 
J sss ^ 


rip 
Jbs 


i 


i (Hi) 


iM 


ff = 


ff = 


ff 



i (iv) I (iii) (iv) | | (iv) (iv) | 

J.0 ji"") r(l> (™) r<p r<p^ (™) rfp'^ 

Jb2b ^ J ssb ^ Jss Jss Jb2 

(EH) I 

J ss 

Several of these form factor relations and consistency checks have been presented in this 
paper. The proof of further relations will be published elsewhere. 
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Appendix 

A Derivation of properties of generalized form fac- 
tors 

In this appendix we derive the formulae for form factors of section |^. We use LSZ tech- 
niques [0 (see e.g. [Q) and assume in addition "maximal analyticity" which means 
that all singularities originate from physical intermediate states. For simplicity we con- 
sider only particles with the same mass m. Generalizations to the case of particles with 
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different masses are obvious. 



A.l Properties of generalized form factors for the pure bosonic 
case 



As usual we write the in-field as 

dp 



{ai^{p)e-^''- + ap{p)e'^''). (A.l) 



27t2uj 

It fulfills the Klein-Gordon equation (9^ + m^) = and when acting on states of the 
form it creates anti-particles and annihilates particles. The commutation rules of 
the creation and annihilation operators are 

[<'ip') , <ip)] = (A.2) 
[<'ip') , <Hp)] = Sa'a2u27rS{p'-p)=Sa>a4n6{d'-e). (A.3) 

Corresponding formulae hold also for the out-field. 

For the matrix elements of a local scalar operator O = O{0) we have the LSZ-reduction 
formulae |57] 

(A.4) 

= % {■■■,p'i I <\pi) o\...Y:.+ij d\ X {...,p[\T [ojI (x)] I ...y: e-'^^- 



(... rf^x \T[Oji,{x)\ \p„...)Z 

where T is the time ordering operator and the source term j{x) = {d"^ + m^) is given 
by the interpolating field (f){x). If pi {p'l) corresponds to an anti-particle (a particle) has 
to be replaced by j. On further reductions and combined with the assumption of maximal 
analyticity the LSZ-formulae imply the crossing formula (|3.4] ) for the connected part of 
the matrix element. We call a contribution to a matrix element °"*( . . . | (9 | pi, . . . )*" 

• disconnected with respect to pi, if its support as a distribution with respect to pi is 
point like and 

• connected with respect to pi if it is a boundary value of an analytic function of the 
Mandelstam variables sij. 

With this notation the first terms in eqs. (|A.4|) are disconnected and the second ones are 
connected with respect to pi or p[, respectively. 

If we interchange in eqs. ([A.4|) in and out the time ordering is replaced by anti-time 
ordering. Comparing eq. (|3.3| ) with (|3^) this means that s + ie is replaced by s — ie. 
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Combined with the completeness of the in- and out-states we obtain the general Watson's 
equations ( |3.8|) . However, from integrability follows a stronger formula. To show this, 
we consider the branch point S12 = {ttii ± 7/12)^ separately. For simplicity we assume 

mi = 7712. 



Lemma A.l Let the S-7natrix factorize as deTioted i7i ( \3. and let s\2 he in a neigh- 
bourhood of Aw? or and all other Sij away from Am? , then 



Fa [si2 + ie, + ie 



F^,{si2-te,s,, + te)S^^^^lis,2) for S12 ^ W (A.5) 



F^{si2 + ie,Sij + ie) = F^(si2 - ^e, + ie) for S12 ^ (A. 6) 



with a = («!, 02, . . . , On) and a' = {a[, a'g, . . . , a„) and 2 < i < j < n. Corresponding 
formulae hold for all other branch points Sij = Am? . 



Proof: By means of formula (^^^ we may cross all particles except 1 and 2 to the left 
hand side. Using again LSZ we have for the full matrix element 



out 



out 



{p3,...,Pn\ 0\pi,P2y 



(A.7) 



out J 



P3, 



,Pn\a 



out t / 



{pi)0\p2)+7 / d^X°^'{ps,...,Pn\T\Oj^ 



\P2)e~ 



■ipix 



{p3,...,Pn\0\pi,P2) 



out 



(A.8) 



out I 



{P3,...,pn\a^''\p,)0\p2)-i i d'x°-'{p3,...,PrAT* Oj\x) \p2)e 



-ipix 



where we have omitted the indices a and T* means anti-time ordering. The term 



out I 



P3, 



Pn I a°"* ' (pi) 0\p2) is disconnected as in eq. ( |A.4[ ), whereas 



°"*( J93, . . . ,p„ I a™^(pi) 0\p2) in general contains also connected contributions. However, 
for factorizing S-matrices this term is given by 

,Pn|a™^(Pi)ki,...,gn.)*"^"(gm,...,gi|0|P2) 



out I 



P3, 



which is disconnected with respect to pi. Therefore, if we take the connected parts of 
eqs. ( |A.7D and ( |A.8| ) we obtain as in equations ( |3.2| ) and 



°"*(P3, ...,Pn\ 0\pi,p2 )Tonn. = F^iSu + ^C, (trs) (l<r<2<s<n.) , {Skl + ^c) (2<fe</<n.) ) (A. 9) 
°^'{P3, ...,Pn\0\pi,P2 )ZL. = F^{Sl2 - ie, {trs){l<r<2<s<n), {Skl + ^e) (2<k<l<n)) (A.IO) 



which implies the first claim. Moreover, crossing in equations ( |A.7| ) and ( [A.8| ) in addition 
also particle 2 to the left hand side, by the same arguments we confirm the second claim, 
since | pi = | pi )°"*. 
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As a consequence of this lemma and the bose statistics of the particles we have the 



property (i) (c.f. eq. (|3.1CI|) ) 

f^.ij...i- --i^i^ • • •) = ■■,0j,0i,.. .) Sij{9i - 6j). 



(A.ll) 



Iterating this formula we find that for 6i < . . . < 6n the auxiliary function yields 
the matrix element for an out-state 



fZ..a^ {Oi,...,en) = {On, • • • , ^i) s<::< (0 = ( 1 o I Pi, . . . , p„ 



out 
a ■ 



(A.12) 



We obtain property (ii) (see eq. ( ^.11| )) by comparing the analytic parts of the crossing 
relations of the following lemma. 

Lemma A. 2 In terms of the auxiliary functions crossing for the full matrix elements 
reads 



l{pi\0\p2,...,Pn)t.n 

f l{pi\p2)2ftn{03,...,er,) 

+ C''fZ.niO^+^7r.,...,9^) 

l{Pl\Pn)nfZ.n-liS2, ■ ■ ■ , On-l) 
f2i'..nl(^2, • 



(A.13) 



+ /f..„i(^2,...,^n,^i-^vr_)Cii 



for 9i>e2>...>e.n 



for ^2 > • • • > On > 9i 



where 7r_ 



TT — e. 



Together with property (i) this lemma implies the general crossing formulae ( |3.14| ) for 
arbitrary ordering of the rapidities. 

Proof: The disconnected contributions in eq. (|A.13|) follow directly from the LSZ - formu- 
lae ([A.4|) . Moreover the LSZ - formulae imply that crossing of particle 1 means pi — > —pi- 
In terms of the Mandelstam variables s or the rapidity differences (see also figure ^ for the 
analytic properties) this means sij + ie tij — ie or 6ij in — 6ij. However, since there 



is no branch cut separating region II and IV (see (|A.6|) ), this is equivalent to considering 
sij + ie ^ tij + ie or 6ij ivr + Oij. Hence, because of 6ij = \9i — 9j\ we have the 
equivalences 

9ij ^ in + 9ij <^ 9i ^ 9i + in for 9i > 9j 
9ij in + 9ij 4^ 9i—>9i-in for 9i < 9j 

which imply the claim. 

The form factors have poles determined by one-particle states in any subchannel 



a 



J, . . . , LXj 



a,) C (ai. 



if the square of the total momentum of all particles in the 
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subchannel equals the one-particle mass squared. We follow the arguments of []T3[ and 
in particular of (see also [0, ^). A particular type of poles is always present, even 
if there are no boundstates. These poles are often refered to as kinematic poles. If for 
instance in ( |A.13D particle 1 is the anti-particle of 2, then one-particle intermediate states 
with the quantum numbers of all other particles j (2 < j < n) yield contributions to a 
pole at pi ~ P2, since then {p2 + Pj — Pi)^ ~ rri'j in eq. ( [A.13| ), (see figure 0). The residue 



Figure 7: A graph contributing to a pole of a form factor. The dashed lines denotes an 
off-shell line with the propagator -, — : ^-^^2 — 2 ■ 

of this pole is given by property (iii) of eq. ( ^.12 ). This can be seen as follows. By (ii) we 
have for 61 ^ 62 and 61 6j, (j = 3, . . . , n) 

Oi-ie- O2 

/f..„l(...,01-m_) ^ ] Ci2^?(g2,---,gn) 

Oi + ie- 02 

for some function (7(6*2, . . . , 6*^). Employing now the well known identity = ^^i7i6{a), 
with P denoting the principal value, the general crossing relations ( |3.14| ) imply for the 
full matrix element for 61 ^ 62 and 61 ^ 6j, (j = 3, . . . , n)) 

l{pi\0\p2,...,Pn)Z.n ^ l{Pl\P2)2fZ.n+[j^^+^^m-02))c,2g 

^ l{Pl\P2)2 fl.n S2n ■■■82,+ [j^^ " ^^'^(^l " ^2)) C12 g. 

Comparing the delta-function parts we obtain 

^(^2, . . . , ^n) = 2^ /3^..„(^3, . . . (1 - 52n ■ ■ ■ ^23) 

which yields (iii). 

If there are also bound states there are additional poles and we have also property 
(iv). The latter is obvious from figure |^ and eq. ( p.lO[ ) up to a normalization. The 
normalization follows from the following argument developed in [^]. Let us consider a 
model with a bound state of type c of two particles of type a and h such that the attractive 
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region is connected analytically (by a coupling constant) to a repulsive region, where the 
bound state decays. For simplicity we consider first the two-particle form factor 

fZ{ea,) = {0\O{0)\pa,Pb)2, iOa>e,) (A.14) 

such that the (scalar hermitian) operator O connects the bound state c with the vacuum 

/f (^e) = ( I O{0) I = ^ 0. (A.15) 

Then in the atractive region of the coupling the two-point Wightman function reads 

( I 0{x)0{y) I ) = Z^A+(x - y, ml) + contributions from other masses (A.16) 

or the time ordered two-point function in momentum space fulfills 

(0\ d(k)O(0)\0) ^ Z'^— atk^f^mj (A.17) 

-ml + ze 



where is a wave function renormalization function. In the repulsive region the contri- 
bution from the two-particle intermediate states ab is given by 



{o\o{x)o{y)\o) = /^^^(o|0(^)|p-p^>)ir"(p-p^>|0(y)ior'+ 

= ^j_jefS,{e)f''\e)A^{x-y,s,,) + ... (a.is) 

where summation over the multiplets a and b is assumed with 5^5 = m^+m^+2mamft cosh 6 
and A_|_(x, m^) = (27r)~^/ cP'ke^^^Q{kQ)2Tx5{k'^ —m?). In the repulsive region the functions 
fah{(^) aiid f'^"'^{9) have poles in the unphysical sheet at i^^^ (Im 9'^^^ < 0), respectively. 
If we move to the attractive region these poles will cross the integration path and by 
analytic continuation we get 

(0|O(x)C»(y)|0) = ^(-/ +/ + r|rf^/»r>) A4x-y, .„,) + ... 



ab ab 



^ [Res - ^Res ) f^,{0)f'\e)A4x -y,ml) + ... (A.19) 



ab ab ) 



Both residues give the same contribution, because 

Res a{0)r\B) = Res f^^{-e) S^,{9)r\e) = fg{-9l^v':R,^l,r\ei^ 

ab ab 

Res fZ{e)r\e) = Res /„^(-^)5.„(-^)r'"^(-^) = -/„^(-^SJ</:'?^c<^L/'''"(^:J 
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where property (i) eq. ( p.lO| ) and the residue formulae for the S-Matrix ( |2.10| ) and ( |2.11|) 
have been used. Using this and comparing eqs. ( |A.1(: ) and (|A.19 ) we obtain 



which agrees with (iv). The general case may be proven similarly. 

A. 2 Properties of form factors for the general case 

We now consider the case in which the particles are taken to be fermions and the operators 
may be of fermionic or bosonic nature. Again we use LSZ techniques and "maximal 
analyticity" . The two component fermionic in-field is 

= / ^ {<iP) <P) + v{p) e^^^') (A.20) 

and fulfills the Dirac equation {i'yd — m) il){x) = 0. The anti-commutation relations are 

{a:,%p'),a^\{p)} = 

Corresponding formulae hold for the out-field. We use the conventions for the 7-matrices 
and the spinors of ( p.29| ) and ( |5.30| ). 

The LSZ-reduction formulas for fermions read 



out 



l{...,p[\o\p,,...)z... 

+ tld^x :^^{...,p[\T[Oj^,ix)] I ...)^>(j9i)e-^Pi^ (A.22) 
: ao.; °!*(... \Oa:)^{p[)\p^,...)Z... 

-zaoa'Jd'x \t\Oj^,^{x)] \p,, . . .)Z..v{Pi) e'"'^' (A.23) 



where j(x) = ip{x) {i'yd + rn) and aoa = — 1 if C is fermionic and apa = 1 otherwise. 
Obviously aoa = ("1)" if n is the total number of fermions in the states. A similar 
formula holds if we interchange particles and anti-particles. The invariant form factors 
G'*^'-'^(sjj+ie) defined by eq. (|3.13|) are again boundary values of analytic functions. Again, 



if we interchange in and out time ordering is replaced by anti-time ordering, which means 
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again that Sij + ie is replaced by Sij — ie. The crossing relation for the connected part of 
the matrix element reads 



out 



m 

{-ir n ^Oa. E ^(Pn) ■ ■ ■ U{pm) T u{pm+l) " " " v{p,) . . . pf^ (A.24) 

O 



i=l 



Watson's equations for the invariant form factor functions G acquire the same form as 



those for F for the bosonic case (|3.^ ). Also Lemma |A.1| holds for the invariant form 
factors G. 

Analogously to the bosonic case it is convenient to introduce the vector valued auxiliary 
function f^iS.) which is considered as an analytic function of the rapidities of the particles. 
Its components coincide again with the physical matrix elements for a particular order of 
the rapidities. 



f^ie^,...,e^):={o\o\pu...,p^y: , for e,>...>er^ 



(A.25) 



In the other sectors the function f^{6) is again given by analytic continuation. Again, 



as a consequence of Lemma A.l and the fermi statistics of the particles we have now the 



property (i) in the form (c.f. eq. ( p. 17] )) 



f...ij...{- ■ • • • •) = f...ji...{- ■ • • • •) {-S)ij{Oi - 6j). 



(A.26) 



The LSZ-formulae ([A.22|) and (|A.23| ) imply the general crossing formulae ( p.l8| ). Note 
that some signs in these formulae depend on the choice of the relative phases of the u- 
and the f-spinors taken in eq. ( p. 301 ). The crossing formulae again, as for the bosonic 
case, implies the properties (ii) and (iii) as given by ( |3.17 ), where we have used the fact 
that aoi = o"2n ■ ■ ■ <^23, if particle 1 has the same statistics as particle 2. 



B Proof of the theorem [^7l| 



In the proof of theorem ^]T| we follow ||5^ (see also ||53|| ) 



Proof that eq. fulfills (i) 



Property (i) (c.f. eq. (|3.17|) ) follows directly from the Yang-Baxter equations, the defini- 
tions of the soliton-soliton S-matrix ( p.l5| ) and the pseudo-ground state fl and Watson's 
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equations for F{9) 

FiOji) ^...ji... C,,ji„X- ■ ■ ^jj di- ■ ■) Sij{6ij) = F{6ji) SijiOij) C,„ij...{.- ■ ■ di, ■ ■ •) 

= —F{9ji) a{9ij) C„,ij...{. . . 9i, 9j . . .) 

= F{9ij) ^...ij... C,„ij„X- ■ ■ 9j . . .). 

The minus sign is due to the fermi statistics of the sohtons (c.f. (|0|)) 



Proof that eq. ( gj^ fulfills (ii): 

Using (i) the property (ii) (c.f. eq. ( |3.17| )) may be rewritten as a difference equation 

fZ.n{0) = fZ.nmQl...n{0) (B-l) 



where 9^ = {9i, . . . ,9'„ = 9n — 2tii) and Q{9) is the trace of the monodromy matrix (^4.1| 
over the auxiliary space for the specific value of the spectral parameter 6*0 = 9n 



Qi..M) = ^^oTIM) , with Tl,M = Ti..M,0n) 

since S'„o(0) = P„o is the permutation operator. This may be depicted as 



(B.2) 




In the following we will suppress the indices 1 . . .n. The Yang-Baxter relations 
the well known commutation rules for the matrices A, C, D defined in eq. 



imply 



0(9, u)C{9, v) 
C{9,u)A{9,9) 

C{9,u)D{9,9) 



C{9, v)C{9, u 
a{9-u) 



h{9 - u] 
aiu — 9] 



A{9, 9)C{9, u) - ^^^—^A{9, u)C{9, 9) (B.3) 



F>{9,9)C{9,u) 



b{9 - u) 
c(u - 9' 



D{9,u)C{9,9) 



' b{u-9) ' ' b{u-9y 
In addition there are commutation rules where also the matrices A^, C^, defined by 



1 [^u.) - y^Q^Q^^ D«(0j 



are involved [O 



C{ff,u)A'^{9_) 
C{9:,u)D'^{9) 



''\AQ{9)Ci9,u) - 5^ ''l Ai9:,u)CQ{9) (B.4) 



b{9n - u) 



u 



b{9n - u 
c{u - 9'J 



b{u 



j{D{&,u)C^{9) 
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To analyze the right hand side of eq. ( [B.l| ) we proceed as follows: We apply the trace of 
T'^d) to the CO- vector /'^(^') as given by eq. ( [4.10| ) and the Bethe ansatz ( [4 .41 ). In the 
contribution from A^{6_) 



nc{e;,ui)---c{e!,ur 



e. 



Ul 



Ur 



because of charge conservation only the amplitudes a{9'^ — uj) appear in the S-matrices 
S{6'^ — Uj) which are constituents of the C-operators. Therefore we may shift all Uj- 
integration contours Cq' to Ce_ without changing the values of the integrals, because the 
functions a{6'^ — Uj)(t>{9'^ — uj) are holomorphic inside Cq' — Cg. 

We now proceed as usual in the algebraic Bethe ansatz and push the A'^i^O) and D^{ff) 
through all the C-operators using the commutation rules (|B.4|) and obtain 



U ^7?^^ A'^iS.) C{9, Ul) ■ ■ ■ C{9, um) + uwA . 

b[en - Uj) 



(B.5) 



■C{e_,Um)+UWD. (B.6) 



fA Ku, - 

The "wanted terms" written explicitly originate from the first term in the commutations 
rules ( [B.4|) ; all other contributions yield the so-called "unwanted terms". If we insert these 
equations into the representation ( [4.10|) of /(^') we find that the desired contribution from 
A^ already gives the result we are looking for. The wanted contribution from applied 
to Q gives zero. The unwanted contributions cancel after integration over the Uj. All 
these three facts can be seen as follows. We have 



QA'^ie) 



i=l 



nD'^{9) = 



which follow from eq. 

The relations ( [4.15 ) for (j){u) and ( 4.121 ) for F{6) imply that the wanted term from A'^ 
yields f{0). The commutation relations ( |B.3| ) and ( |B.4| ) imply that the unwanted terms 
are proportional to a product of C-operators, where exactly one C{6_,Uj) is replaced 
by C^{6_). Because of the commutativity of the Cs it is sufficient to consider only the 
unwanted terms for j = m which are denoted by uw^{u) and uw^{u). They come from the 
second term in (|B.4[) when A^{6_) is commuted with C{6_, Um)- Then the resulting A{9_, Um) 
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pushed through the other Cs and taking only the first terms in ( [B.3| ) into account and 
correspondingly for D^{6_] Um) we arrive at 

- n to) c(...o . . .c«,., 
-s(M) - - to) n c(g'. «.)... c«(i,. 

Using again ( [4 .61 ) and the relations ( ^.151 ) and ( [4.16| ) for and t{u) we obtain 

where also c{u)/h{u) = —c{—u)/b{—u) has been used and v/ = . . . ,u'^ = + 27rz). 
Therefore after integration of the A-unwanted term along and the D-unwanted term 
along Cqi both cancel. 



Proof that eq. ( ^J^ fulfills (iii): 



We will prove that eq. (|4.1CI|) fulfills (iii) (see eq. (^.17|)) in the form of 



Res ^{9i, . . . , On) — —Con 2« Ci„ ^-1(^2, ■ ■ ■ , ^n-l) (l2...n-l ~ " " " 5'„_i„ 

which is equivalent to eq. ( p.l7| ) due to (i). We consider the n-particle form factor function 
given by eq. ( [4.10| ) 



Ul) ■ ■ ■Ci„J9, Ur, 



with the scalar function 



n-1 

i=2 ' ' i=l 

m—1 

X 



n (0(^1 - 0(^n - r{u, - uj) e±^'(2.™-e,-e„) 

where = 02,..., ^n-i and u = Ui,..., Um-i. We calculate the residue of this function at 
di = On + iTT. It consists of three terms 

Res. f^,J0)=R^ + R2 + Rs 

This is because each of the m integration contours will be "pinched" at three points: 
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(1) Uj = 6n = Oi- in, 

(2) Uj = 6n + in = 6i 

(3) Uj = On — in = 6i — 2in. 

Due to symmetry it is sufficient to determine the contribution from the M^-integration 
and multiply the result by m. 

The contribution of (1) is given by Wm-integration along the small circle around u^ = 
On (see figure^). The S-matrix S{On — Um) yields the permutation operator 5(0) = P and 
^(^i — Urn) the annihilation-creation operator S{in) = K 

7 7 



Caf3 C^'^ — 5s^ 



a 



a 



Therefore we have for Um = On = Oi — in 



^l...n Ci„JO, Ui) ■ ■ ■ Ci,„n(0, On) 



Oil 



u. 



m "n 



Oo 



On-1 



Ui 



Ur. 



Or, 



m—1 



= n {H^^l - Uj)a{On - Uj)^ Cin ^2... n-lC2...n-l{0.,Ul) ■■ ■ 

..n-l{0_, Um^i) S2n ' ' ' Sn-ln 

where Ci„ is the charge conjugation matrix with Cq/j = dap- We have used the fact 
that because of charge conservation the amplitude h{-) only contributes to the S-matrices 
^(^^i — Uj) and a(-) to the S-matrices S{On — Uj). 

We combine this with the scalar function and after having performed the remaining 
Mj-integrations we obtain 



R\ 



(1) 



i:7&~ n %s {-2m) Res a(6'„ - Um)4>{0n - Um) h{Oi - u„)0(6'i - m„) 

iV„ 9 dl=e„+fK Um=dn 



*n-2 

n-1 



X 



F{in) n {f{0^ - 0,)F{0, - On)m 



— Ur 



i=2 

m—1 



— Uj)T[Uj — U„ 



m—fl—On) 



2i aOnCin f2...n-l{Q) " " " Sn-ln 
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if we relate the normalization constants by the recursion relation 

,(/r"(o))' 



n 



n-2 

-1)" and 



47rm 



(B.7) 



We have also used that aon 

F{u) F{u + in) (f){u) = 1 and b{u + in) (f){u + in) a{u) (f){u) r{—u) = 1 

which follows from the definitions ( [4.12|) and ([4.13|) . Finally we have used the normal- 
ization F{in) = 1 and ReSu=Qa{u)(f){u) = ReSu=inb{u)(f){u) = — 2z//™"(0) because of 
eq. (|4.14| ). Note also that the signs from the as and 6s cancel and aon = since 
all particles are fermions. 

The remaining contribution to (iii) is due to R2 and 



p(2) p(3) 
-n-l...n "T ^l...n 



2« Cln /^..n-l 



If both particles at 1 and n are solitons both vanish. If one particle at 1 or n is an soliton 
and the other an anti-soliton one term gives the desired expression and the other vanish. 
If both particles at 1 and n are anti-solitons both terms cancel. These fact can be proven 
as follows. 

The contribution of (2) is given by the Wm-integration along the small circle around 
um = (^1 (see again figure ^). Now S{9i — Um) yields the permutation operator S'(O) = P 
and the co-vector part of this contribution for Um = 61 = 6n + in is 

fil...n Ci,„n{d, Ui)--- Ci,„n{0, Um = Oi) P„{s) 











Ui 








Ui 


































01 




V 






Um 


















9^ ^02 










n 
i=l 




m—l 

-Um) n 




- Uj)d{6n - Uj)) 








X 


Cln ^2...n-lC2...n- 




lui) ■ ■ ■C2...n-liO.,U 


m-l) Pl{s) 


Pn 


(s) 





(B.8) 



where the Yang-Baxter relation (p.8|) has been used iteratively. Pi{s) and P„,(s) project 
onto the components where the particle at 1 is a soliton and at n is an anti-soliton, 
respectively. The remaining components if both particles at 1 and n are anti-solitons are 
calculated below. 
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We combine this with the scalar function gn and after having performed the remaining 
Mj-integrations we obtain 

>(2) 



xm 



Res (— 27ri) Res d{6i — Ur, 



n-1 

xF{in) n (^(^1 - ^i)Fi9^ - en)a[ti, - u 



1 - Mm ) aWn - Ur, 



i=2 



m—1 



n(- 



±s(2«„-6»i-9„) 



= -2^ao„Ci„/2^..„_i(£)Pi(s)P„(s) . 

if we apply the condition exp(27rs) = (—1)" = (Ton and if we relate the normalization 
constants as above. We have used the identities 



F{-u)F{u + i7r)d{u) 



1 and a(u) 



[uj a(u — m I 0(u — ITT I T(—ui 



1. (B.9) 



The contribution of (3) is given by Wm-integration along the small circle around Um = 
On — iiT (see again figure |^). Now S{9n — u.m) yields the annihilation-creation operator 
S{m) = K and the co-vector part of this contribution for Um = Oi — 27ri = On — in is 



^l...nCi,„n{9,Ui) ■ ■ Um = On- in) Pi{s) 











Ul 
















Um 


Um 


Oi 


02 


^n-l 




\0n 



Oi 



Ur 



'n-1 



J 



Ul 



Ur 



m—1 



— Ur 



n [b{0,~U,)b{On-U,)^ 
i=l 3=1 

X Cin ^2...n-lC2...n-l{0, Ui) ■ ■ ■ C2..,„-l(£, Mm-l) Pl{s) Pn{s) 



(B.IO) 



if the particle at 1 is a soliton and that at n an anti-soliton. This contribution obviously 
vanishes if n is a soliton. Again the remaining components if both particles at 1 and n 
are anti-solitons are calculated below. 

Again we combine this with the scalar function gn and after having performed the 
remaining Wj-integrations we obtain 

Rf„nPl{s) = C,nfZ.n-MPl{s)Pn{s) 
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T7§~ n %S {-2tIi) Res 6(6'„ - M„)0(6'„ - Um) b{9i - - Mm) 



n-1 
1=2 

m—l 



provided that we fix the normahzation constants as above. We have used the identities 

b{u) F{—u) F{u + iTr) (f){u) = 1 and a{u) (f){u) a{u — in) (f){u — in) t{—u) = 1. (B.ll) 

Finally we calculate R^^^n + ^i^ln the case that both particles at 1 and at n are 
anti-solitons. Instead of eq. (|B.8|) we have now for Um = Oi = 9n + in 



Y[d{6i-Um) n [a{6i-Uj)d 



u 



a{9n — Um~l] 



i=l j=l 

X ^2...n-lC2...n-liO,Ui) ■ ' " C2...n-1 Mm-2) ^2...n-l Mm-l) ^1 (s) P„,(s) + . . . 



and because of the Yang-Baxter relations 

Cilui) ■ ■ ■C{lum^2) Dilum^l) 

= n u ~ ^ ^(^' «i) ■ • • «m-2) + . . . 

where the dots refer to similar terms with D2,..n-i{0,Uj), (j < m — 1). Because of 
symmetry with respect to the C-operators it is sufficient to consider only this term. 
Similarly we get instead of eq. ( p.lOp Um = Oi — 2ni = On — in 



W b{9i - Um) n (^(^1 ^ ^j)^ 



i=l j=l 



b{9i - Um-l] 



X VL2...n~lC2...n-l{d,Ui) ■ ■ ■ C2...n-l{0, Um-2) ^2...n-l (£ ^^m-l) Pl{s) P„(s) + 

and because of the Yang-Baxter relations (p3.6| ) 

C{lu,) ■ ■ ■C{lum-2) A{lum-l) 

= n l( I u^^i) cii m) • ■ ■ c{i um-2) +... 
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where the dots again refer to similar terms with A2...n-i{0.,Uj), {j < fn — We apply 
D2...n~i{^,Um-i) and A2...n-i{0,Um-i) to the pseudo-vacuum, use as above the identities 
( [B.9| ) and ( [B.ll| ), and find that the sum i?^^^ + R^^^ is proportional to (with u = Um~i) 

f , I c(9n - u) ■ ,^ , "^^^aiui-u) , 

4* n m - u)m - u) n - «) 

0(t^l - «j j=l ^ %) 

Due to crossing we have 

C{9n - Um-l) _ C{9i - Um-l " 27ri) 

a(6'„ - Um-l) K^i - Um-l - 27ri) ' 

In addition we use the quasi-periodicity properties ( |4.15|) and ( |4.16| ) of (f>{-) and r(-) and 
get 



lc-g+2ni JCg^j h{9i-Um-l) 

TT a{9i - Um-i)4>{9i - Um-l) n TT^ T^K' ~ ^m-l) 

t=2 j=l b{Um-l - Uj) 

which vanishes since the integrand is holomorphic inside the contour Cg — {Cg + 27rz). 

C Some useful formulae 

In this appendix we provide some explicit formulae (which partly may also be found 
elsewhere in the literature) for typical scattering matrices, minimal form factors and 
some auxiliary functions which we frequently employed in the explicit computations. We 
state some typical integral representation, which are very useful since via ( 4.19|) and ( |4.20|) 



they relate the scattering matrix and the minimal form factors effortlessly. The infinite 
product representations in terms of Gamma functions, obtained from the evaluation of the 
integrals or the direct solution of the functional relations, make the singularity structure 
more transparent. For numerical purposes it is often more useful to express the Gamma 
functions with the help of Euler's product representation in terms of rational functions 
at the cost of an additional infinite product. 
A typical S-matrix eigenvalue is (for a > 0) 

a + X dt 



bvnix,a) = = exp / — 2e smhra; 

a — X Jo t 
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According to ( |4.19| ) and ( |4.20| ) the corresponding minimal form factor function is therefore 
f {z7ix,a) = exp / — 2e 



t 2 sinh t 

{2l + 2 + a-x){2l + a + x) + f ) 



E {21 + 1 + ay ' r(i + f-f)r(f + !)• 

In particular, for a = we recover the scattering matrix of the Ising model 

S = -\ /™"(i7rx) = sin^x 

For negative values of a we use S{Q,a) = 1/S{9, -a) and f'^'''{9,a) = -a). A 

further typical S-matrix eigenvalue^ is (for < a < 1) 

^, , sinf(a + a;) dt sinht(l-a) 

S nnx, a) = 4- = exp / — 2 smh tx . 

sm ■2(0 — 0;) JO t smht 

with the corresponding minimal form factor function 

rminr N /-"^ sinht(l - a) 1 - cosht(l -x) 

/™"-(z7rx,a) = exp / — 2 



nn 



t sinh t 2 sinh t 

2l + 2k + a + x 2l + 2k + 2 + a-x f2l + 2k + 3-a 



k=oi=o 2l + 2k + 2-a + x2l + 2k + A-a-x \2l + 2k + l + a^ 



n 



r(A; + 1 - f + f ) r(A; + 2 - f - f ) (T{k + | + f 



2 



.=0 r(fc + f + f) r(A; + i + f-f) vr(fc + |-f; 

The sine-Gordon soliton-soliton S-matrix reads 

dt sinhi(l-z/)t 



bssU'n'x) = annx) = exp / — — j— smhte 

JO t smh i^ut cosh i^t 

00 00 

nn 



2 ^^'j*^ 2 

2/ + z/ + fcz/ + x 2/-h2 + fci/ + x 2l + l + iy + kiy-x 21 + 1 + ku - x 



^,^qi^q21 + V + kv - X 21 + 2 + kv ~ X 21 + I + V + kv + X 21 + I + kv + X 
^ T{{v + kv- x)/2) r((2 + kv~ x)/2) r((l + v + kv + x)/2) r((l + A:;/ + x)/2) 
Ik + fci^ + a;)/2) r((2 + kv + x)/2) r((l + i/ + fcz/ - x)/2) r((l + fci/ - x)/2) 

with 

TT 7r 

a(z7r(a; + z/)) = — cot —x cot —{x + v) a^inx). 
Consequently, the minimal form factor function is 



■^For instance almost all diagonal scattering matrices related to perturbation of certain conformal field 
theories may be built out of these elementary blocks. 
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/--(^TTx) = exp r 
Jo 



°° dt sinh 1(1 1 - coshi(l - x) 



2 

t sinh ^vt cosh |i 2 sinh f 



OC' OO OC' 



-j-j- j-j -j-j- 2m + 2Z + 2 + i/ + A;i/ — x 2m + 2Z + i/ + /ci/ + x 



fe=o /=o m=o + 2Z + 3 + i/ + /ci/ — X 2m + 2Z + l + i/ + /ci/ + x 
2m + 2Z + 4 + ^i/ - X 2m + 2Z + 2 + + X 



2m + 2Z + 3 + A;i/ - x 2m + 2Z + 1 + + x 

/2m + 2/ + 2 + z/ + A;z/ 2m + 2/ + 2 + A;z/\ ^ 
\2m + 2l + l + p + kv 2m + 2l + 3 + kuJ 

^ ^t{1 + ^{3 + u + ku - x)) T{l + ^{l + iy + ku + x)) 

ioio r(Z + 1 + + - x)) r[l + + + x)) 

r{l + l{3 + ku-x)) r(/ + i(l + A;i/ + , 



r(/ + 2 + i(A:z/-x)) r(/ + l + i(A;z/ + x)) 
r2(/ + l(l + i/ + H) T^{l + l{3 + kiy)] 



r^{i + i + ^{u + ku)) r^{i + i + ^ku) 

with asymptotic behavior for \^9\ oo, {\^6 — tt\ < |(3 + u — |1 — u])) 

(ivr - ^) = (e^l^l + o(l)) 

with the constant 

sinh 1(1 — 1 — i/\ 



1 r"^ dt ( si 

2 Jo T V sinh I 



^vt cosh sinht i/t 



The corresponding functions 0(w) = {F{u)F{m -\-u)^ and T(ii) = (^(f){u)(f){—u)^ with 
F{i'Kx) — sin(|x)/J^*"(i7rx) are 



sinh 1(1 — ( cosht(| — x) — 1 



i?2^Hj sin(7ra;) Jo t sinh | i/t sinh t 

1 ^^2l + l + iy + ku + x2l + 2 + u + ku-x( 21 + ku + \ 



nn 



3 



F'^(^-^^ t=oi=o 2l + ku + x 2l + l + ku-x \2l + u + ku + ^ 

1 ^ r{l{ku + x)) r{lii + ku-x)) r^{^ + i) 



with 



F2(f) ^^\r[l{l + iy + kiy + x)) T[l{2 + iy + kiy-x)) r^{^ + l) 



^ " cosf(x + i/)^^ ' 
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, F\i7^/2)F\-i7,/2) . 

TyiTTX) — ; SIWKX SlUyTTX / V) . 
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